Universal spectral form factor for chaotic dynamics 



Stefan Heusler^, Sebastian Miiller^, Petr Braun^-'^, and Fritz Haake"'^ 
^ Fachhereich Physik, Umversitdt Duisburg-Essen, 45117 Essen, Germany 
^Institute of Physics, Saint-Petersburg University, Saint-Petersburg 198504, Russia 



m 
o 
o 

in 



Q 
U 



> 

(N 
O 
On 
O 

m 
o 



> 
X 



We consider the semiclassical limit of the spectral form factor K{t) of fully chaotic dynamics. 
Starting from the Gutzwiller type double sum over classical periodic orbits we set out to recover 
the universal behavior predicted by random-matrix theory, both for dynamics with and without 
time reversal invariance. For times smaller than half the Heisenberg time Th (x h~l'^^, we extend 
the previously known r-expansion to include the cubic term. Beyond confirming random-matrix 
behavior of individual spectra, the virtue of that extension is that the "diagrammatic rules" come 
in sight which determine the families of orbit pairs responsible for all orders of the r-expansion. 



Introduction: One of the fascinating quantum signa- 
tures of chaos is universal behavior of the correlation 
functions of the spectral density of energy levels, for gen- 
eral hyperbolic dynamics [1]. Three universality classes 
were suggested by Dyson and Wigner; one, called "uni- 
tary" , has no time reversal symmetry, while the other 
two do have Hamiltonians H commuting with an anti- 
unitary time reversal operator T; if = 1 one speaks 
of the "orthogonal" class while the "symplectic" case has 
= — 1. The Fourier transform of the two-point corre- 
lator of the level density, called spectral form factor, is 
predicted by random-matrix theory (RMT) [2] as 



A'u„i(T) = T , K^rthir) =2t-t ln(l + 2r) . 



(1) 



in the unitary and the orthogonal case; here r is a 
time made dimensionless by referral to the Heisenberg 
time Th oc h~-^^^, with / the number of freedoms; 
the results (1) hold in the (semiclassical) limit of large 
dimension of the matrix representation of H and for 
times up to the Heisenberg time, < r < 1. The 
orthogonal form factor allows for the Taylor expansion 
KoMt) =2t- 2t2 2t3 . . ., for < r < i. 

Understanding the observed fidelity of individual dy- 
namics to RMT has been an elusive goal, in spite of 
considerable efforts based on parametric level dynamics, 
semiclassical periodic-orbit theory, and the so-called non- 
linear sigma model [2]. We shall here take a non-trivial 
step towards that goal, on semiclassical ground. 

Periodic-orbit theory a la Gutzwiller [3,2] gives the 
form factor of an individual spectrum as the double sum 
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where S^,T^, and are the classical action (including 
the Maslov phase), period, and stability amplitude of the 
7th orbit; the angular brackets demand averages over (i) 
the center energy {E -\- E')/2 in the product p{E)p{E') 
of two level densities before doing the Fourier transform 
w.r.t. the energy difference E — E' and (ii) over a time 
interval small compared to the Heisenberg time. Most 



orbit pairs 7, 7' interfere destructively in the double sum 
(2). Finite contributions arise only from families of pairs 
wherein the action difference can be continu- 

ously steered through the quantum scale h toward zero, 
by varying parameters defining the family (for early pre- 
monitions see [4]). Periods and stability amplitudes do 
not differ noticeably within such an orbit pair. 

Berry's "diagonal approximation" [5] includes the triv- 
ial pairs {7,7} and, given T invariance, {7,7} where 
the overbar indicates time reversal; it yields if^^-* = 
PT,-y ^) = where /3 = 1 without and /3 = 

2 with T invariance, due to the doubling of contributing 
pairs in the latter case. The sum J^-y l^7p'^('''~ ^) — '''i 
known as the sum rule of Hannay and Ozorio de Almeida 
(HOdA) [6], reflects ergodicity for long periodic orbits. In 
view of (1) the diagonal approximation gives ifuni in full, 
and the first term of the r-expansion of ii'orth- 

Sieber and Richter [7] recently found a family of or- 
bit pairs which for the T invariant Hadamard- Gutzwiller 
model yields the quadratic term of the r-expansion, 
K^l^ = -2r2, as in (1). Each Sieber-Richter (SR) pair 
has a close self-encounter which in configuration space 
looks like a small-angle crossing for one orbit and like a 
narrowly avoided crossing for the partner orbit. Gener- 
alizations to arbitrary hyperbolic systems with two free- 
doms were given in [8-10] and for more freedoms in [11]. 

Before identifying the new families of orbit pairs giving 
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— ^/ , iVjjjjj — we must briefly review SR. 

Each self-encounter involves two orbit stretches which 
are nearly mutually time reversed; it may be depicted 
as ^ (or X ) , arrows indicating sense of traversal. On 
either side of the "encounter graph" =^ , each of the two 
orbits has a long loop attached. Assuming symbolic dy- 
namics available (to uniquely define periodic orbits and 
even, approximately, short orbit stretches by symbol se- 
quences; our results are valid more generally) we could 
write E and E for the two (nearly) mutually time re- 
versed orbit stretches in the encounter region, R, L for 
the two long loops, and L for the time reversed of L; we 
may thus write EREL for an orbit and EREL, LERE 
for its SR partners [12,8]. Note that the orbits in a SR 
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pair traverse one loop in the same sense while the senses 
of traversal are opposite for the other loop; here, T in- 
variance is seen as required for SR pairs to exist. 

Following [13,9,10] we parametrize an cncoimtcr with 
the help of a surface of section V transverse to an orbit, 
say 7 = EREL, somewhere within the encounter; V is 
two dimensional for / = 2, the case we limit ourselves 
to. The stretch E pierces through V in a, point Xa which 
can be made the origin of a coordinate system spanned 
by tangent vectors e., and e„ to the stable and unstable 
manifolds of 7 through Xa in V. A second piercing is 
associated with E and thus opposite in sense; it happens 
after the traversal of the right loop at a point Xb. We 
define a close-encounter region by requiring the unsta- 
ble and stable components of the difference Txb — Xa to 
respect a classically small bound c independent of %, 

Txb — Xa = Ueu + SCs , |u|<c, |s|<c. (3) 

Moving V we leave the encounter after a time t„ given 
(asymptotically) by |u|c^*" = c. Conversely, we find the 
start of the encounter going backwards in time by tg with 
|g|gAts _ jjgj-g ^ jg ^j^g Lyapunov exponent of the sys- 
tem. The duration of the encounter thus is 

^enc = ts +tu = - in . (4) 

A \us\ 

Roughly, linearization of the dynamics about any point 
within the encounter breaks down at either end. 

As was shown in [13,8 10] the partner orbit EREL 
pierces through V first in xP^ = Xa+ uBu and then in x^ 
with Tx^ = Xa + scs- Moreover, the orbits in a SR pair 
differ in action by the area of the parallelogram spanned 
by the four points Xa,xP,Txb,Txl [9,10], 

AS = us; (5) 

that product is canonically invariant and thus indepen- 
dent of the precise location of the surface V (The vectors 
igjCu are pairwise normalized as eg A e„ = l). Inas- 
much as weighty contributions to the form factor must 
have AS = 0{h), we can conclude that the duration of 
relevant encounters has the order of the Ehrenfest time 
Te = jlnj-, much smaller than the period T = 0{Th)- 
To evaluate K'^'^\t) we need the cumulative duration 
P{u, s\T)duds of all orbit stretches within self-encounters 
of a long orbit 7 of period T = 0{Th), with unstable and 
stable components of Txb — Xa in [u, u+du] and [s, s+ds]. 
Ergodicity yields (Refs. [8-10] use different conventions) 

P{u, s\T)duds = T{T - 2Unc)^^^duds (6) 

with D. the volume of the energy shell. This results from 
integrating the ergodic return probability density 
over the two times of piercing of the orbit through a sec- 
tion V. The factor T indicates that one piercing, say the 
one at Xa, may occur at any time in the interval [0,T]. 



The time of the subsequent piercing at xi, can then lie 
only in an interval of length T — 2tenc) hence the sec- 
ond factor in (6); this is because both traversals of the 
encounter region have length fcnc and may not overlap. 
(Overlapping stretches E, E are either impossible, as in 
the Hadamard-Gutzwiller model [12], or indicate an or- 
bit with a self-retracing loop identical with its SR partner 
[8]). Note that in the density P{u,s\T) each encounter 
is weighted with the duration A/'tenc ; the combinatorial 
factor N = 2 arises since in 7 = EREL the two stretches 
E,E are equivalent; we must therefore employ ^^A/t'eiJ^ ' 
to count each encounter only once [14]. 

The contribution to the form factor reads K^^^'^^-^ = 
(E^ l^fl'^r - ^)) j:juds^P{u,.s\T)2cos{us/h). 
The sum ( E^ • • ■ ) gives the factor r through the HOdA 
sum rule, as for K^^\ The integral gets no contribu- 
tion from the leading term T^/fi in P, and this is why 
the 0{Te/Th) correction in (6) is important; the inte- 
gral with P — > —2tenc'P/fl becomes independent of the 
bound c in the limit h ^ and is easily found (since tgnc 
cancels from the integrand and Th = ^/2nh) as — 2r. 
The RMT result /T^^) = _2-,-2 jg thus recovered. 

One might wonder why no "parallel" encounters with 
graphs =C come into play. The simple reason is that 
the would-be SR partner of an orbit EREL decomposes 
into a "pseudo-orbit" (separate periodic orbits ER and 
EL), not admitted to the Gutzwiller sum in the first 
place. However, parallel encounters will be met with 
below. The quantification (3) of "close" must then be 
changed to Xb — Xa = ue-u + se^ , |w| < c , |s| < c. 

The foregoing review of SR has highlighted those twists 
of the original formulation of [7], in particular the ap- 
pearance of ^o'l'^^"^'' ) which make the extension to higher 
orders in r a rather elegant travail, to which we now turn. 

Orbits pairs contributing in third order: We now 
present five families of orbit pairs relevant for t^. They 
can be constructed from two SR "switches" / X . 
Three families have two separate encounters wherein the 
orbits differ and four intervening "loops" near identical 
for both orbits; two more families arise as one of the loops 
shrinks, to let the two encounters overlap, see Fig. 1. 

Starting with two independent encounters, we must 
obviously check three possibilities: (aa) both "antiparal- 
lel" , pictorially =^ =^ , (ap) antiparallell and parallel, 
i.e. ^ =t , and (pp) both parallel, i.e. =£ =t . For 
all of them there are two distinct ways of filling in inter- 
vening lines, with the two encounters either in series (s) 
or intertwined (i). In symbolic notation, those two ways 
read for case (aas) E1AE1BE2CE2D (in series) and for 
case (aai) E1AE2BE1CE2D (intertwined); here Ei,Ei 
{E2,E2) represent the two orbit stretches of the first (sec- 
ond) antiparallel encounter which are nearly mutually 
time reversed, while the remaining symbols refer to in- 
tervening lines; one immediately checks that only the aas 
orbit allows for a partner orbit, E1AE1BE2CE2D, while 
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the would-be partner of the aai orbit turns out a pseudo- 
orbit decomposing into the two periodic orbits E1AE2C 
and E1BE2D. In the same vein we find that for {ap) 
only the api orbit and for {pp) only the ppi orbit have 
non-decomposing partners. We have thus identified three 
families of orbit pairs with two far apart self-encounters. 
Clearly, ppi pairs arise even without T invariance, while 
aas and api pairs do require that symmetry. 




ac. pc 

FIG. 1. The five orbit pairs entering t^. Labels describe 
encounters as antiparallel or antiunitary-symmetry-required, 
parallel, intertwined, serial, and cloverleaf. 



Now on to the families resulting from shrinking away 
one of the four loops in the foregoing families. The three 
remaining loops make for a cloverleaf (c) structure; the 
encounter region accommodates a triple of oriented short 
{0{Te)) stretches with encounter graphs ^ {pc) and 
^ (ac); in symbolic notation, they involve E,E,E and 
E, E, E, respectively. Schematically, the two types of 
cloverleaf orbits look like the thick lines in Fig. lac,pc. 
We shall argue that each has a unique partner, shown as 
dashed lines in Fig. 1. The family (pc) with three paral- 
lel orbit stretches ^ in the encounter does not require 
T invariance while the ^ family (ac) does. 

To check the uniqueness of the partner for ac and pc 
wo start from thc^ respective thick lines in Fig. 1. Each 
candidate for partnership must have its three long loops 
nearly identical with those of the starting orbit, save pos- 
sibly for time reversal; those loops are differently con- 
nected in the encounter region. A sextet of cloverleaf 
orbits thus seems to arise, whose encounter graphs are 
(without arrows, momentarily) = ^ ^ >< >< . 
The last three immediately drop from candidacy since 
they represent SR partners of = contributing to or 
entail decomposing pseudo-orbits. Now putting arrows 
on = , say for ac as ^ and hooking on the three loops 
as in the thick line of Fig. lac, we find only ^ to lead 
to a non-decomposing partner, the dashed line in Fig. lac 
with the encounter graph ^ (plus, of course its time 
reverse); the uniqueness of the pc partner (up to time 
reversal, if T invariance holds) is shown similarly. 

Duration and action differences: For the independent- 
encounter pairs aas, api, ppi the action differences of the 
two encounters sum up to AS' = ?ii.si + U2S2- The triple 
encounters ac,pc require extra thought. 



Beginning with ac we employ a surface of section V 
through the encounter region and denote the three points 
of piercing by Xa,b,c- As in (3) we have Txb — Xa = 
UbCu + SbCg, Txc — Xa = UcCu + ScRg] thc bouud c must be 
respected by all six distances |wf,|, |s&|, \uc\, \sc\, |u6— Wc| = 
\ubc\, \ sb — Sc\ = Sbc and the cloverleaf encounter lasts for 



X max{|u6|, |mc|, |w6c|}niax{|s6|, |sc|, |s6c|}" 

Again, relevant encounters will have durations of the or- 

1 ^ 

der of the Ehrenfest time ts = j\n^. 

Now look at an ac pair {7, 7^*} with encounter graphs 
^ and ^ . To find the action difference we proceed in 
two steps. In the first, we employ an auxiliary orbit, the 
SR partner 7' of 7 related to the encounter of the "upper 
two" stretches, labelled o, b, leaving the "lowest" stretch 
c as in 7. For 7' the encounter region is ^ . Accord- 
ing to what we have said above about the four points of 
piercing through the surface of section of an SR encounter 
the piercings of 7' occur at x'^ = .x,, + u^eu and xj, with 
Tx'j^ = Xa + SbCs] moreover, the action difference within 
{7,7'} is given by (5) as S^' — = UbSb- In our second 
step we arrive at the ac partner 7^ of 7 as thc SR partner 
of 7' w.r.t. the stretches a, c. Once again invoking what 
we know about the four piercings in an SR encounter we 
have x^ = x'^ + [tic — U}y)eu and Tx"^ = x'^ + .ScC.,, and 
the action difference S^p — S^i = {uc — Ub)sc- The action 
difference of the ac pair (valid also for pc) thus reads 

A/S = S-yP — Sj = UbSb + UcSc — UbSc ■ (8) 

The term UbSc represents an "interaction". 

-contributions from two simple encounters: We first 
generahze the density (6) to P{ui, Si, U2, S2IT), (up to the 
factor duidsidu2ds2) the "area" of times (ia^^a) ^ [0,?"]^ 
such that the points {x^ = x{t^), ^ = 1,2} of piercing 
are followed by piercings at .x^ = x{t'^) with unstable 
and stable components of x'^ — (for parallel encoun- 
ters) and Tx'^ — x^ (for antiparallel encounters) in the 
intervals [u^, + du^], [s^, + ds^]. We obtain 

P{U1,SI,U2,S2\T) = ^T(T -2(^encl + ^enc2))^ (9) 

by integrating the ergodic probability density fl^^ for two 
encounters over the four times {^ai^b j A* = 2}, respect- 
ing the order of those times dictated by the "diagrams" 
aas, api, ppi and the general rule that an orbit must leave 
one encounter region before reentering or entering the 
next one; the latter rule in fact separates independent- 
encounter from cloverleaf families. The restrictions on 
the times of piercing in question give rise to the small 
but decisive corrections icnc,^ in (9), as before in (6). 

Again in analogy with (6), the density (9) overcounts 
a pair of simple encounters, by a factor which we shall 
argue to be A/'tenciienc2- Obviously, the times of pierc- 
ing may lie anywhere during the respective encounters, 
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hence the product of the two durations. The factor M 
is of combinatorial nature. For instance, in aas pairs 
the two antiparallel encounters are indistinguishable such 
that Naas = 2; likewise, J^api = 2 since the two stretches 
of, say, the parallel encounter are indistinguishable; fi- 
nally, in ppi pairs all four orbit stretches (symbolically, 
El, E2, El, E2) are indistinguishable, hence Af^ 



For T invariance, (2) gives 



(3) 

orth,ind 



ppi 



= T J_^dui . . .dS2 



— ^P('Ui . . . S2|T")2cos-^ where we have already 
allowed the HOdA sum rule to yield a factor t as above 
and used the overcounting factor Af as an indicator for 
the cases aas,api,ppi. Using the action difference (5) 
for both encounters and the weight (9) it is easily found 
that only the part gjp- x 3 x 4 x T'^tenci^enc2 in P sur- 
vives the integration, and actually yields the square of 
the twofold integral met with for r^; all other terms and 
all c dependence vanish with h ^ 0. We thus have 



^(3) ^ ^ 3 ^ 

orthjind ji^ 



_ j 4t'^ for aas, api 
\ 2t^ for ppi . 



(10) 



For dynamics without T invariance, however, only ppi 
pairs exist and yield K^^^ = r^, one half the ppi term 
in (10) since a ppi orbit now has no time reverse. 

-contributions from triple encounters: Both for ac 
and pc pairs, by reasoning as before we have the density 



P^\ub,Sb,Uc,Sc\T) = 



T{T-Mli,+...)' _ 3TH 



enc 



2^2 
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where the dots point to terms killed by the integration 
to come, as ?i ^ 0. The crucial term oc tf^^ is due to 
minimal loop lengths. The orbit must leave an encounter 
before reentering in the antiparallel or parallel sense; oth- 
erwise, an SR pair already accounted for in t"^ would arise 
in the antiparallel case, whereas the parallel case would 
lead to a new family with stretches involved in encoun- 
ters resembling multiple repetitions of shorter periodic 
orbits; such families turn out irrelevant for U ^ 0. 

To count each cloverleaf only once we divide out the 
familiar Mf^^^. For pc encounters we have Afpc = 3 since 
the three parallel stretches 3 are indistingushable while 
ac encounters ^ entail Afac = 1- Given T invariance, 



the form factor picks up K 



(3) _ 
orih.cl 



P'=\ub . . . Sc\T)2cos^. The limit h 



. dut, 
yields 



-^orth.cl 



AT' 



~6t^ for ac 
-2t'^ for pc . 



Without T invariance, no ac pairs exist and pc pairs, not 
accompanied by time inverses, give but K^-^ ^.j = — r^. 

Conclusions and outlook: Adding contributions from 
independent and cloverleaf encounters we get K^,.^i^ = 
(4-1-4-^2-6 - 2)r3 = 2t^ and K^J- = 0, as in (1). To 
third order in r at least, then, semiclassical treatment of 
individual hyperbolic dynamics gives the universal form 
factor characteristic of ensembles of random matrices. 



The five families of orbit pairs met here resemble di- 
agrams known from field theoretic treatments of disor- 
dered systems [15] and from quantum graphs [16]. The 
analogy between classical orbits and diagrams in field 
theory should persist in higher orders. Orbit pairs (alias 
diagrams) with n—1 separate simple encounters con- 
tribute to r" ; upon shrinking intervening loops we expect 
to find all other relevant orbit pairs. The weight of each 
family includes a correction, due to the ban of encounter 
overlap and small as a power of ^ , exclusively affecting 
the form factor. 
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